Effective Hamiltonian for fermions in an optical lattice across Feshbach resonance 
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We derive the Hamiltonian for cold fermionic atoms in an optical lattice across a broad Feshbach 
resonance, taking into account of both multiband occupations and neighboring-site collisions. Under 
typical configurations, the resulting Hamiltonian can be dramatically simplified to an effective single- 
band model, which describes a new type of resonance between the local dressed molecules and the 
valence bond states of fermionic atoms at neighboring sites. On different sides of such a resonance, 
the effective Hamiltonian is reduced to either a t-J model for the fermionic atoms or an XXZ model 
for the dressed molecules. The parameters in these models are experimentally tunable in the full 
range, which allows for observation of various phase transitions. 

PACS numbers: 03.75.Fi, 67.40.-w, 32.80.Pj, 39.25+k 



Recently, there are many exciting advances in the ul- 
tracold atoms physics For these developments, 

two experimental control techniques play the critical role: 
one is the Feshbach resonance to control the interaction 
magnitude between the atoms and the other is the 
optical lattice to introduce diverse interaction configura- 
tions 01 . It is natural to consider combination of these 
two techniques, and indeed, significant efforts have been 
put forward towards this direction |E S III ISl ISl ■ 

A fundamental problem along this direction is to de- 
rive an appropriate Hamiltonian for this strongly inter- 
acting system which can serve as the starting point for 
further investigations. A number of generalizations of 
the Hubbard model have been proposed to describe this 
system by including the on-site atom-molecule coupling, 
typically ignoring the upper-band occupations ^ (see 
the comment in Ref. 7]). These generalizations may 
model physics associated with a very narrow resonance, 
however, they are not adequate to describe typical broad 
Feshbach resonance, such as for '^^K or ^Li. For the latter 
case, first, one needs to include all the multiband coupling 
terms even if the system temperature is well below the 
band gap. The reason is that the strong atom-molecule 
coupling results in population of the upper bands (the 
on-site coupling rate is typically larger than the band 
gap) Second, one also needs to include the atom- 
molecule coupling from the neighboring sites which has 
been ignored in all the previous works. We will see that 
off-site atom-molecule coupling is typically larger than 
the atom tunnelling rate, and inclusion of these off-site 
interactions leads to qualitatively different physics. 

In this paper we rigorously derive the interaction 
Hamiltonian for fermionic atoms in an optical lattice 
across a broad Feshbach resonance, taking into account of 
both multiband couplings and off-site interactions. The 
strong on-site interaction between the atoms make them 
first form local dressed molecules. Under typical experi- 
mental conditions, we then derive an effective single-band 
Hamiltonian, describing the resonant interaction between 
the local dressed molecules and the valence bonds (sin- 
glets) of fermions at neighboring sites ^3- this ef- 
fective single-band resonance model (different from the 
conventional Feshbach resonance), multi-band occupa- 



tions have been incorporated through the local dressed 
molecule states. On different sides of this resonance, the 
effective Hamiltonian is reduced to either a t-J model for 
the fermionic atoms or an XXZ (anisotropic Heisenberg) 
model for the dressed molecules, opening up the prospect 
of using this system to probe some fundamental physics 
associated with the latter two models. 

We consider fermionic atoms with two internal states 
labelled by the spin index cr =t,i. The atoms are 
loaded into an optical lattice, and tuned close to a 
Feshbach resonance by an external magnetic field. 
The Hamiltonian, in terms of the field operators 
vj/(ni) ^ xj/t^") respectively for the bare molecules and 
the fermionic atoms, then has the form H = Hq + Hj, 
with Ho = E.^T,iI'^i"^HTa + Va)^i^U'r + 
/*(")t(T„^ + y„^+j,j^)^Md3r and Hi 

In the above expression, the kinetic energy Ta — 2Tm — 
—h^y^ /2m (m is the atom mass), and the potential 
energy Va = = Vq [sin^ kox + sin^ koy + sin^ koz] . 

The Va and Kn are due to the optical lattice potential 
from far-off-resonant laser beams with a wave vector 
ko = 27r/A. For simplicity of the notation, the potential 
depth Vq is assumed to be the same along the x,y,z 
directions. The detuning vi, of the bare molecules can 
be controlled by an external magnetic field B. The 
atom-molecule coupling rate a and the atom back- 
ground scattering rate Ubg are determined from the 

atom scattering length as a = \J Anh^^coW \ab\ /m, 

Ubg = 4:Trh?ab/'rn, where we have assumed the atom 
scattering length near Feshbach resonance takes the form 
as — ab{l — W/ {B — -Bo)), with Oh, the background 
scattering length, Bq, the resonance point, W, the 
resonance width, and /ico, the difference of the atom 
magnetic moments between the closed and the open 
scattering channels. 

The filed operators ^'i"'' and can be expanded 

with the Wannier functions associated with the lattice 
potential in the forms 5'^"'' (r) = 



binWo"' (r - r,), where Wp ' (r - r^) 
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wfj {x - Xi) wfj {y - yi) wfj {z - z,) (w^' (r - r^)) 
are the Wannier functions for atoms (molecules) at the 
site with p = (PxiPyTPz) labelling different lattice 
bands, and ajpo- (bip) are the associated mode operators. 
With these expansions, the Hamiltonian Hq then has the 

form Ho = J2^p (4'"'' + f^b) 4pb^p + E^p<x ^P^alpa<^'P<y + 

E^p EjeAr(i) (t'^^blpb^p + ^ alpa<^jP<^) ^p to 

the nearest-neighbor tunnelling, where N (i) de- 
notes the neighboring sites of i, — e'^'' w 
[2{px + Py + Pz) + ^ ^/VqE^ under the harmonic approx- 
imation to the potential well {Er = h?kQ/2m is the atom 
recoil energy), and the tunnelhng rates t'^\tp^ can be 
determined through the standard band calculation (see 
Fig. 1 and Ref. [l|). 

In the expansion for Hi, usually one ignores all 
the terms except the ones for the on-site interaction. 
However, for typical broad Feshbach resonance, the 
nearest-neighbor atom-molecule coupling rates cj"™^, = 
a J Wp™^* (r) wi"'^ (r) ui^"' (r — 6r) (Pr {\Sr\ = \rj — ri\ is 



the lattice constant) can be significantly larger than the 
atom tunnelling rates tp"^ , and thus should not be ne- 
glected. To see that, we have numerically calculated 
the exact Wannier functions and their overlaps, and 
some of the results are shown in Fig. 1. The calcu- 
lation clearly shows that the neighboring couplings are 
not negligible. For instance, for the lowest band with 
p = s = s'= (0, 0, 0), the rates cj"™^, and tp ^ scale with 
the potential depth Vo/Er by roughly the same expo- 
nential form, and the ratios between them are estimated 
to be c[''^p^Jt'^^ « IQVWe^ {lA^/^JE'r) respec- 
tively for ^Li (^"if) (see Fig. 1). For this estimation, 
we have taken the following parameters for ^Li {'^'^K): 
W = 180G (8G), at = -2000aB (ITOas), A = l^m 
(0.8/xm), ^co ~ {^J-B■ Bohr magneton; as- Bohr 

radius). So, we include here all the nearest neighbor cou- 
pling terms in the expansion for Hi. The next nearest 
neighbor couplings can be safely neglected unless the lat- 
tice potential is very weak with Vo/Er < 3. Up to this 
order, the interaction Hamiltonian Hi is expressed as 



Hr = 



ipqs 
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[%■,ps^'^pa^s^a,^'l + h 
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ipqss' jeN{i) 



pss' pqss' 



■i;pqss'a»pia!qT [o-is^ajs'i " ai^iajs']) + C2;pss'Oipa3^] a.js'i + h.^ 



pss' 



(1) 



where Cn^^L = " / Wp™''* (r) wi"^ (r) w^f^ (r) d^r, 4°^s' = " / w^^* (r -I- 5r) w^"^ (r) w^?^ (r) d^r, 4°pqss/ = 



''Oipss' 



f/,,/<^*(r)<)*(r)z.(^)(r)u;lf)(r)d3r, c^^;^^, ee U,, J wl^> (r) w^^^* (r) wi^^^r) wif^r ~ Sr) d\ cg^,,, ^ 
Ubg J Wp^* (r) Wq^* (r) wi"'^ (r - Sr) w^"'' (r - Sr) d^r, 4?pqss' = i^bg) J Wp^* (r) Wq^* (r - 6r) wi""^ (r) w^"-' (r - Sr) (fi 



The above Hamiltonian seems extremely complicated, 
however, we show now that under typical experimen- 
tal conditions it can be dramatically simplified. To 
make this possible, first we note there are several 
different energy scales in this system, including the 



(am) 
'^0;pss' 1 



the band 



2^/VoEr, the off-site interaction energy 



on-site interaction energy E, 
gap Ehg ^ 

^off ~ , and the atom tunnelling rate Et 

In experiments, the energy scale set by Eon,Ebg is 
typically much higher than the one set by Eoff,Et. So 
we take the approximation {Eon, E^g) 3> {Eoff,Et). In 
this case, we can first solve the single-site Hamiltonian 



H, 



Ep [(4' 



blpb^P 



E, 



(a) t 



Epqs (4?^UlpaisTa»s'i + /i-c.) 



Epqss' 4'?^qss'°ipi"IqT°^«Ta*s'i on the site We 
assume that the average atom filling number n of the 



lattice is smaller than 2. If a site has a single atom 
on the pth band, its energy is just given by ep"^-*. If 
two atoms occupy the same site, they will form a local 
dressed molecule state, which in general can be written 
in the form l^&i^j) — |0), where |0) denotes the 
vacuum state and the dressed molecule creation operator 
4p = Ep X^cpblp + Epq iMPqoIpxalqT • In this expression 
for d|^, the superposition coefficients X/ipj7Mpq' with 
the normalization J^p X^pX^'p + Epq 7;^pq7/.'pq = 5^,^,,, 
are determined by sol ving the Schrodinger equation 
Hi l^f^^) = E^ \^if_,) 01, where fi labels different 
eigenstates with the corresponding eigen-energy £^^. 
This kind of two-particle equation has been solved 
in Ref. with a Harmonic approximation to the 
potential well. Here, we only need to mention two 
general features of the local dressed molecule states: 
first, the eigen-energies of the dressed molecule can 
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FIG. 1: The normalized atom tunnelling rate tn = t'^'' /Er 
(the subscript means the lowest band), and the normalized 
overlaps between the atomic and the molecular Wannier func- 
tions uq = -y/ A/2 J ^Wq"''' {x)j ^'"'o"'' (a;)^ dx (on-site) and 

ui = / (w,)™' (a;)) * 4"' (x) {x + A/2) dx (neigh- 

boring sites) shown as a function of the lattice depth Vo in 
the unit of the atom recoil energy Er- The "o","x 
points denote the exact results from the numerical calcu- 
lation, and the dashed curves show the fits from the for- 
mula t„ ~ (3.5/0F) {Vo/Erf^^ exp (-2y'Vo/Er^ , uq ~ 

0.77(Vo/£^)'''^ ui ~0.52 {Vo/Erf^'^exp (^-2y/Vo/Er 



be tuned by the external magnetic field B through their 
dependence on i/^, although such a dependence is in 
general nonlinear; second, the energy difference between 
the adjacent eigenvalues is typically of the order of 
the band gap energy E^g. 

We consider the case with one of the eigen-energies, 
say Ef_ig, pretty close to the two-atom free energy 2epa 
on a certain band po, i.e., we tune the magnetic field 



B to satisfy the condition 



< E, 



bg. 



illustrated in Fig. 2. The and po can be chosen re- 
spectively as the ground state of the dressed molecule 
and the lattice lowest band, although they are not sub- 
ject to such a restriction. Under the above condition, 
if the atoms start in the band po, their state evolution 
will be restricted in the Hilbert subspace involving only 
the excitations of the modes d],, and a]^ „, as all the 
other states are significantly detuned by a energy scale 
of {Eon,Ebg) which is much larger than {Eoff,Et) 
So each site can only take four possible states, given by 
|0), \a) EE al^, |0) {a =T, i), and \d) ee d^^^ |0). We can 
then project the full Hamiltonian H = Hq + Hi into the 
physical subspace specified by the projector P = ^ Pi, 
withP.^lO), (0| + |T),(T| + |i),ai + M),(d|. After such a 
projection, the effective Hamiltonian i?e// = PHP takes 
the form 



Hr 



ff 



^A{B)4d, + Y, ^ 



i jeN(i) 



E E 



^(^gdl (aiT^ji ~ aa^it) + ^-c-) + tdPd\djP 

ta ^ Po,l^O,jaP + tda ^ C?" dja^-^a^^ + CdTlidrijd + Ca {niTlj/A- • Sj) 



(2) 



where we have used the simplified notations di = d. 



A (B) EE E^„ (B) - 2e^°^ (the atom 



energy zero point has been shifted to coincide with epj), and defined the number and the spin operators 



riid 



- J 



d'-di, Ui 



J2cr°'icrO'i'j, Si = Z^crcr' "L^'^'^'"iff'/2 wlth thc PauH matrix a = {a^, ay, az)- The param- 



eters in the Hamihonian H^ff are given by g = 7^popo*PoV2 +J2qX 

X)q IX^ioql +2 Re(^qgg, X^j^qC2.qsg/7;^os's) + X^pqss' 7^iopq'^2;pqss'7/Jns's, 



(am) 
Moq''l;qPo 



Po X/qs 7, 



qs /Moqs l:qspoPo' 



X/q ITmo 



qpo 



td = 

|2 ,(a) 



-2Re(V X 



qs A;^oq'^l:qP( 



(aa) 
'^3;poPoPoPo ■ 



,qs7a'ospo) 2 Re(^pqg 7*|^pqC 



Jaa) 



llpqPosT/Jospo <^d — Z^pp'qq'ss' ')'AioPP''''Moq'q''3;pqss''7posp'7/Joq's' I 
I 



The Hamiltonian H^ff represents an effective single- 
band model, but it has incorporated multi-band infor- 
mation into its parameters through the structure coef- 
ficients XMopi7Mopq of the local dressed molecules. This 
Hamiltonian contains no on-site interactions as they have 
been exactly taken into account into the structure of the 
dressed molecules. As one tunes A (B) through the mag- 
netic field, the Hamiltonian H^f / describes a crossing res- 
onance between the local dressed molecules di and the 
atomic valence bonds ai^aji — aiiaj^ on the neighbor- 



ing sites. This resonance may have a richer structure 
than the conventional Feshbach resonance between the 
bare molecules and the free atoms as first, the dressed 
molecule here has been a complicated superposition of 
the bare molecules and the local Cooper pairs, and sec- 
ond, the resonance to valence bonds on neighboring sites 
may introduce rich configurations depending on the lat- 
tice geometry. 

The Hamiltonian i?e// supports rich physics. Its de- 
tailed study will be presented elsewhere. Here, we investi- 
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FIG. 2: Illustration of resonance between the local dressed 
molecules and the atom valence bonds on neighboring sites. 



gate Heff in two limiting cases with the detuning |A {B)\ 
significantly larger than the neighboring coupling rate g. 
We can see already rich phase diagrams from -ffe// in 
these limiting cases. First, we consider the case with 
the population dominantly in the atoms and the local 
dressed molecules only virtually excited due to the large 
detuning A (B) (the lattice filling number n < 1 in this 
case). We can then reduce i?e// to an effective Hamil- 
tonian involving only the atomic operators ai^. For this 
purpose, we define the projection operator Pa = P/"' 
withP/"^ = |0),(0| + |t),(T| + U),;ai, and keep the terms 
up to the order of \g\^ /A{B) in the projection. The 
reduced Hamiltonian can be derived through either the 
second-order perturbation or a canonical transformation 
[ll| . and it takes the form of the famous t-J model 



tj 



E 

i;jeN{i 



ta Pa0.iryO,jaPa 



J {Si 



run-j/A) 



H: 



(3) 

where the parameter J = — Cq — 2\gf /A{B) In 
Htj^ we have not included the 3-site nonlinear tunnelling 

{diiaji — flijaji), which is usually omitted in the liter- 
ature 17]. The t-J model plays a fundamental role in 
study of high Tc superconductivity 0, 0|. De- 
pending on the ratio ta/J, the atom filling n, and the 
lattice dimension and geometry, the t-J model shows 
rich phase diagrams, including, for instance, the atomic 
anti-ferromagnetic phase, the d-wave superconductivity 
^M, and likely the resonating valence bond (RVB) states 
|1(I ITll |. The t-J model was previously derived as the 
strong interaction limit of the Hubbard model |0, ^| . 
Here, although our basic Hamiltonian H^ff does not re- 
semble any forms of the Hubbard model, we get exactly 
the same reduced Hamiltonian as the t-J model; further- 
more, our atomic realization of the t-J model is not sub- 
ject to the constraint J ^ ta as is the case from the 



Hubbard model. All the parameters here can be exper- 
imentally tuned in the full range, which, together with 
the easy control of the lattice dimension and geometry, 
makes this system the ideal test bed for many outstand- 
ing predictions and hypotheses associated with the t-J 
model. 

We now consider another limiting case of our basic 
Hamiltonian Pe// with the population dominantly in 
the dressed molecules and the atoms only virtually ex- 
cited (still, |A(P)| ^ \g\). In this case, we define 

the molecule projection operator P^ = ■ p.^™^ with 

2 



j(m) _ 



|0), (0| + \d)^ {d\. Up to the order of \g\^ /A (P), 
the reduced Hamiltonian for the dressed molecules has 
the form 



H„ 



E 



A'n,rf ■ 



E 



tdPmd\djP, 



(4) 

where t[i = t^ + 2 |.g|' /A (P) , c;, = - 2 /A (P), 
and A' = A (P) + 2nc \g\^ /A (P) [tic is the lattice 
coordination number). This Hamiltonian describes 
hard-core bosons with neighboring interactions, and 
it is identical to the magnetic XXZ model Hxxz 

(id/4) [BeffZ, + E,eiv(.) (X.X, + Y,Y, + S,Z^,Z,) 
{S, = cyt'^, Beff =2{A' -\-ncc'j) /t'^) for the effective 
Pauli operators Xi,Yi, Zi through the following map- 
ping: X, = P„ (^d, -\- dfj P„i, Yi EE iP„ (d^ - dfj P^, 

Zi = 2d\d, - 1 |i,[li|. The phase diagram of the XXZ 
model is know pretty well. For instance, it has the fer- 
romagnetic, the canted XY, and the anti-ferromagnetic 
phases ,19], which correspond respectively to the Mott 
state, the superfiuidity state, and the checkerboard 
state (one occupation every other sites) of the dressed 
molecules. It also has more exotic quantum phases 
(RVB spin liquids etc.) if the latti ce g eometry induces 
frustration of the above spin orders [23. 

In summary, we have derived the effective Hamiltonian 
for fermionic atoms in an optical lattice across broad Fes- 
hbach resonance, taking into account of both the multi- 
band configurations and the direct neighboring interac- 
tions. In certain limits, this Hamiltonian is reduced to 
either the t-J model for the fermionic atoms or the XXZ 
model for the local dressed molecules. The latter two 
models connect to many fundamental physics associated 
with many-body systems. 
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